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behavior in one space dimension exactly. For superohmic dissipation, wo find 
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scaling. The dynamical critical exponent depends on the spectral density of the 
dissipative baths. We also discuss the Griffiths singularities, and we determine 
observables. 
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1. Introduction 

Real physical systems are seldom bereft of quenched disorder that occurs in the guise of 
impurities, defects or other lattice imperfections. This disorder can have a significant 
impact on the physics near phase transitions. These effects are particularly dramatic 
near zero-temperature quantum phase transitions (QPTs) where they give rise to 
a bevy of fascinating phenomena such as quantum Griffiths effects [1] [2], infinite- 
randomness critical points featuring activated dynamical scaling [3l H], as well as 
smeared phase transitions [5] (for recent reviews, see, e.g., [6j[7]). 

It is by now well understood that the enhanced disorder effects at quantum 
phase transitions can be attributed to the disorder being perfectly correlated along 
the imaginary time axis which becomes infinitely extended at zero temperature. A 
particularly important role is played by the rare regions, droplets that are caught 
in the wrong phase due to atypical disorder fluctuations (for instance, ferromagnetic 
droplets in a bulk paramagnet). Due to the perfect disorder correlations in imaginary 
time, these rare regions are infinitely extended in the time direction which severely 
hampers their dynamics and thus leads to the exotic phenomena mentioned above. 
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A dissipative environment further enhances these disorder effects. An example 
is provided by the magnetic quantum phase transition in the random transverse- 
field Ising chain. In the absence of dissipation, this system features an infinite- 
randomness quantum critical point separating the ferromagnetic and paramagnetic 
ground state phases. It displays unconventional activated dynamical scaling 
and is accompanied by power-law quantum Griffiths singularities [Sj [4]. Ohmic 
dissipation prevents sufficiently large ferromagnetic droplets from tunneling [8l [9l 
[To] . This destroys the sharp quantum phase transition by smearing, as was first 
predicted by phenomenological arguments [S] and later confirmed by strong-disorder 
renormalization group calculations [11] [121 US] ■ 

Recently, some of us (M] [15] considered the influence of Ohmic dissipation on the 
quantum phase transition of a disordered rotor model with continuous 0{N) order 
parameter symmetry. Such a model appears as the order parameter field theory of 
the itinerant antiferromagnetic transition [16] or the superconductor-metal transition 
in nanowires |17| . By applying a strong-disorder renormalization group, it was shown 
that the critical behavior of this transition is controlled by an infinite-randomness 
critical point that falls into the same universality class as the random transverse-field 
Ising model. While this result may appear surprising at the first glance (because the 
two systems have different symmetries) , it follows from the fact that the rare regions in 
both systems are exactly at the lower critical dimension of the problem, in agreement 
with the general classification scheme put forward in Ref . |18l . 

The case of Ohmic dissipation and the resulting exotic infinite-randomness critical 
point should be contrasted with the dissipationless case, relevant, e.g., for dirty bosons 
[19j . which features a Kosterlitz-Thouless phase transition with conventional (power- 
law) dynamical scaling. The qualitative difference between these two phase-transition 
scenarios immediately poses the question: What happens in the presence of dissipation 
that is weaker than Ohmic and thus "interpolates" between the two cases mentioned 
above? 

In this paper we therefore study the quantum phase transition of a disordered 
0{N) rotor model in the presence of super- Ohmic dissipation. We apply a strong- 
disorder renormalization group to the order-parameter field theory of the problem. In 
one space dimension, this theory can be solved analytically, giving the asymptotically 
exact critical behavior. Our paper is organized on the following grounds: In Sec. [21 we 
introduce the order-parameter field theory and discuss the phase diagram. Sec. [3] is 
dedicated to the derivation of the strong-disorder renormalization group fiow equations 
which are then analyzed and solved in Sec. 31 In Sec. [SI we compute the behavior of 
a few key observables. Finally, we conclude in Sec. [6] 

2. Order-parameter field theory 

Following Landau, one of the general techniques for studying quantum phase 
transitions consists in deriving an effective field theory in terms of the order parameter 
fluctuations from the the appropriate microscopic Hamiltonian. For quantum phase 
transitions in Fermi liquids, this approach was pioneered by Hertz and Millis [TBI 120] . 
Our starting point is the resulting quantum Landau-Ginzburg- Wilson free energy 
functional (or action). 




(1) 
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Here, </? represents an A^-component vector order-parameter in d space dimensions. 
X = (x, r) stands for both the position vector x and imaginary time r with the 

integration measure J dx = J d'^x /g^^^ dr. u is the standard quartic coefficient. In 
the absence of disorder, the two-point vertex r(x, y) can be defined via its Fourier 
transform which reads 

r(q,w„) =r + ^oV+7ol^n|'^'°, (2) 

where r represents the bare distance to criticality, and is a microscopic length scale. 
The non- analytic dependence of r(q, cj„) on the Matsubara frequency lOn is caused by 
the coupling of the order parameter field to a dissipative bath. We parameterize the 
dissipation by the mean- field dynamical exponent zq. If the dissipation is caused by 
simple oscillator baths of spectral density £{oj) ~ w'*, one obtains s — 2/ zq. We are 
mostly interested in the case of superohmic dissipation corresponding to 1 < zq < 2. 
However, for comparison, we also consider the dissipationless case, zq = 1 and the 
Ohmic case, zq = 2. (Even weaker damping, formally corresponding to zq < 1, is 
irrelevant because the resulting frequency term is subleading to the regular term 
arising in the Landau expansion). The coefficient 70 derives from the the coupling 
of the order parameter to the dissipative baths. We now consider quenched disorder 
which does not change the functional form of the Landau-Ginzburg- Wilson functional 
([ij qualitatively. Instead, the distance from criticality r becomes a random function 
of the position vector x. Disorder also appears in the coefficients ^^0,70, and u. 

Landau-Ginzburg- Wilson theories of the type defined in ([T]) appear in several 
important physical systems. For instance, the case of an 0(3) order parameter 
with Ohmic damping in d = 3 dimensions represents a quantum anti-ferromagnetic 
transition in itinerant electron systems [16| , while an 0(2) order parameter in one 
space dimension is relevant for the superconductor- metal transition in nanowires |17) . 

In the following, we focus on the case of one space dimension and study the 
critical behavior of the model as a function of the dissipation exponent zq. We start 
by analyzing the phase diagram. It is well known that in sufficiently low dimensions, 
the effects of fluctuations are so large that they destroy any long-range order, as was 
rigorously shown by Mermin and Wagner |21| . To determine for which values of zq 
our model can display long-range order, we calculate the strength of the transverse 
fluctuations in a putative ordered phase. Provided that the dissipation remains active 
in the ordered phas^l, the strength of the transverse fluctuations can be estimated by 

where A is a frequency cutoff. In the undamped case, zo = 1, the integral is 
logarithmically divergent in the infrared limit. This implies that true long-range order 
is impossible but suggests a state with quasi-long-range order. In contrast, for zq > 1, 
the integral is infrared convergent, implying that a long-range ordered phase can exist. 

Having established the phase diagram, we now turn to the critical behavior of 
the order parameter field theory ^ in the presence of quenched disorder. We aim 
at applying a strong-disorder renormalization group scheme. As this is a real-space 
based method, we discretize the continuum action ^ in space while imaginary time 
remains continuous. Following Ref. [15], we define discrete spatial coordinates Xj 
and rotor variables ipji^r) that represent the average order parameter in a volume V 

I This is a nontrivial point. In some cases, the action JlJ is only valid in the disordered phase, while 
a gap opens in the spectrum of the bath modes in the ordered phase. 
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large compared to the microscopic scale but small compared to the true correlation 
length ^. In addition, we first consider the large- TV limit of the action ([TJ; we will later 
show that the results remain valid for all continuous-symmetry cases. The discretized 
version of the model in the large- limit and in 1 + 1 dimensions reads as 

5 = ^ (r,; + A, + 7^ k„|'/'°) m^n)? 

i Un 

-T^^(f)i{-UJn)Ji(t)i+l{uJn) ■ (4) 

Here, ri, > and the nearest-neighbor couplings Jij > are random, and 
(j)j{uJn) — Jq^"^ (pj{T)e'^'^"''^dT is the Fourier transform of the rotor variable. The 
Lagrange multipliers Xi enforce the large- constraints (t))^) = 1 for each 

(k) ■ 

of the order parameter components ipi at each site i; they have to be determined 
self-consistently. The renormalized local distance from criticality at site i is given by 
€i = ri + Xi. Note that in the disordered phase, all > 0. 



3. Strong Disorder Renormalization Group 

We now implement a strong-disorder renormalization group for the discrete large- TV 
action (|4]). This method is based on successively integrating out local high-energy 
degrees of freedom [3l IH [22j [231 El]- The initial steps (until the derivation of the 
recursion relations) closely follow Ref. [15], we will therefore only outline the salient 
points of the derivation. 



3.1. Single- cluster solution 

A single independent rotor (or cluster) can be described by the action 



5ci = rE 



(r + A + 7|c^„p/"") |0(c^„)r (5) 



The Lagrange multiplier A is determined self-consistently by the length constraint 
which at zero temperature reads 

1 r°° 1 

1 = - ^ / , (6) 

where e = r -I- A is the renormalized distance from criticality. In the super-Ohmic 
case, zq < 2, the integral can be carried out straight-forwardly, and solving for the 
renormalized distance from criticality yields 

e ^ ^^o/(^o-2) ^ 

This power-law relation between e and 7 (which is proportional to the cluster size) 
shows that a single cluster is below the lower critical "dimension" of the problem. In 
contrast, a high-frequency cutoff A is needed to carry out the constraint integral in 
the Ohmic case zq — 2. The resulting dependence of the gap on the damping constant 
7 is exponential, e = jAe~'^"', signifying that a single Ohmic cluster is marginal, i.e., 
right at the lower critical dimension of the problem. 
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3.2. Recursion Relations 

In the action the and Ji are the competing local energy scales; they 

are independent random variables with distributions 7r(J) and /o(e). In each 
renormalization step, we identify the largest of all of these energies, — max(ei, Ji), 
and then decimate the corresponding high-energy mode perturbatively. This 
procedure relies on the disorder distributions being broad and becomes exact in the 
limit of infinitely broad distributions. For now, we assume that our distributions are 
sufficiently broad; we verify this condition a posteriori. 

Suppose the largest of all energies is a local distance from criticality ("gap"), 
say e2, then we integrate out the rotor if 2 in a second-order perturbation expansion 
with the un-perturbed action given by Sq = 2^Si^„(^2+72 |w„|^^^'')|(/)2(wn)P, and the 
perturbation given by = -TY,^J[Ji4>i{~^n)(t>2{'^n) + J2(t>2{~^^n)(t>?,{'^n)]- This 
calculation is identical to that of the Ohmic case [15] because the leading terms do 
not depend on the value of zq. It leads to the term S = —T^^ J(/3i(a;„)(/93(— a;„) in 
the renormalized action coupling sites 1 and 3. The renormalized interaction strength 
J is given by: 

.1='^. (8) 
£2 

Note that this recursion relation is independent of the exponent zo which controls 
the strength of the dissipation. Its multiplicative structure is fixed by second-order 
perturbation theory. 

Suppose now that the dominant energy in the system is an interaction, say J2 
coupling the sites 2 and 3. In this limit, the two rotors ip2, and </33 are essentially 
locked together and can be replaced by a composite rotor variable ("cluster") 1^2- The 
resultant effective distance from criticality 62 for this rotor can be obtained by solving 
the two-site problem involving the rotor variable ip2 and (^3, 

^(r, + A,+7,|c^„|'/"°)|0,(w„)|2-r^J202(-c^„)03(wn) • (9) 

W„ 1=2,3 I^n 

exactly whilst treating the coupling to the rest of the sites as perturbations. This 
calculation is straight forward but lengthy, details will be shown in the appendix of 
this paper. In the superohmic case, zg < 2, the resulting recursion relation for 62 
reads: 

_ f^Y ' 
J2J V2J \e3, 



const (10) 



Here, the quantity x = (2 — zo)/zo vanishes in the Ohmic limit and takes the value 
1 in the dissipationless case. The extra constant approaches —1 in the Ohmic limit. 
The new damping parameter is simply the sum of the two old ones, 

72 = 72 + 73 • (11) 

The additive structure of the recursion (jlOp in the superohmic case needs to be 
contrasted with the multiplicative form e = 2e2e3/ ^2 arising in the Ohmic case [HlllSj . 
It is also worth mentioning that in the dissipationless case, zq = I, the recursion (|10p 
agrees with that found in the dirty boson problem [19], as expected. 
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3.3. Renormalization group flow equations 

The renormalization group steps outlined above are now repeated many times, 
thereby gradually reducing the maximum energy il. As a result, the probability 
distributions 7r( J) and p(e) change with f2. The renormalization group flow equations 
for these distributions can be easily derived following Fisher's treatment of the random 
transverse- field Ising model 4 . This leads to 

^^^^^ - [irn - Pn) A J) + Pn I dJ, dJ^ 7r( Ji) 7r( J2) 5 (j - ^\ , (12) 



— ^ " ~ ^^^^ + ""JI / dei de2 p{ei) p{e2) S (e""' - ej^'^ - e^"" - const)(13) 

Here, ttq and pn denote the values of the distributions at their upper cutoff fl. The 
flow equation for 7r( J) is identical to the Ohmic case. However, the flow equation for 
p{e) is modified to reflect the additive recursion relation (jlOp . 

We now perform a set of variable transformations for i7, J and e 



fin n 1 

il J X 



1 



(14) 



such that the recursion relations depicted ^ and ([TU)) become simple additions. Note 
that the transformation for (3 turns into a logarithm in the Ohmic limit x — >■ 0. Recast 
in the new variables, the flow equations for the distributions P(C) and i?(/3) read 

dP{0 



dT 



^ + iPo- Bo) P + Roj dCi dC2 P(Ci) P{C2)S (C - Ci - C2) , 



(15) 



dR(B) , .dR , 

- (1 + - (^0 -Ro-x)R 



+ Po f dpi dl32 R{pi) R{I32)6 (fi-pi~p2~ 9(x)) (16) 



where g{x) is a constant of order one. Pq and i?o denote the values of the distributions 
at zero argument. In the Ohmic limit, the flow equations ([T5|) and (|T6l) reduce to those 
derived in Refs. [HI [15] while in the dissipationless case {x = 1) they become identical 
to those of Altman et al. [19] with g(\) = 1. 

4. Solutions of the flow equations 

The complete set of fixed points of the flow equations P3|) and ([TB)) can be found 
utilizing the methods developed in Refs. [HdS]. We have performed such an analysis, 
details will be published elsewhere [55] ■ Here, we instead use the following ansatz [iSl 
for the functional forms of the distributions P(C), and R{(3): 

P(C)=Po(r)e-^«(r)^, P(/3) =Po(r)e-««(r)'' (17) 

This exponential ansatz fulfils the fiow equations (ITSI) and (IT51) for the distributions 
P(C) and R{P). The problem is thus reduced to flow equations for the parameters 
Po(r) and Po(r) which read 

^=-RoPo, (18) 

^ = Po(x-Po). (19) 



Influence of superohmic dissipation on a disordered quantum critical point 7 




Figure 1. Schematic of the renormalization group flow in the Pq-Rq plane. For 
C < Cc, the flow asymptotically approaches a line of stable fixed points that 
represent the ordered phase. In the disordered phase, C > Cc, the flow is towards 
Po 0, and _Ro — > oo. For C = Cc, one flows into the critical flxed point at 
Pq = X, and i?Q = 0. 



It is interesting to note that these equations take the famous Kosterhtz-Thoulcss form 
[27] for all x > 0, i.e., for the entire range of superohmic dissipation strengths. (In 
fact, by a simple rescaling, the flow for any x > can be mapped onto the x ~ I 
case.) A schematic picture of the resulting renormalization group flow in the Pq-Rq 
plane in shown in Fig. [T] 

There is a line of fixed points with Rq — and Pq arbitrary. These fixed points 
are stable for Pq > a; and unstable for Pg* < x. The renormalization group trajectories 
in the Pq-Rq plane can be found by combining (fT8|) and (|T9| . This yields 

Ro = Pq-xIuPq + C. (20) 

The integration constant C controls in what phase the system is. For values of C that 
are larger than a critical value Cc = —x + lux the flow is always towards Pq — >■ 0, and 
Rq — ?> oo. This implies that we asymptotically flow into a regime with J <^ e, i.e., 
the disordered phase. In the opposite regime, i.e., when C < Cc, the flow is towards 
a line of stable fixed points with Pq = finite, and Rq = 0. This corresponds to the 
ordered phase with e J. The critical flow, corresponding to C = Cc, separates the 
two regimes. It terminates at the critical fixed point located at Pq ~ x, and Rq = 0. 

We now linearize ([TSl) and (|19p around the critical fixed point, to study the 
renormalization group fiow in more detail. The critical flow (C = Cc) then takes the 
form 

i?o(r) = ^, Po(r)-x + ^. (21) 

The density of surviving clusters at scale F can be found by integrating the loss due 
to decimations along the renormalization group flow. This yields 

A^(F)=iVoe-^'o'^"'[^°^''^+^"^"'^] (22) 

By substituting the solutions (1^ into this equation, one finds that the number of 
surviving clusters decreases as iV(r2) ~ with decreasing cutoff energy ft. The 
typical distance £ between these clusters thus increases as £(!)) ~ il"^. This implies 
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that the dynamical scahng at the critical point is of power-law type. The dynamical 
exponent z = l/x interpolates continuously between the dissipationless case, z = 1, 
and the Ohmic case whose activated dynamical scaling formally corresponds to 2: = oo. 

It is important to realize that close to the critical fixed point, the interactions J 
are much larger than the gaps e. Therefore, almost all renormalization group steps will 
merge adjacent clusters, while no clusters are eliminated. The number of sites (original 
rotors) contained in the surviving clusters therefore remains essentially constant with 
r 00. As a result, the typical size (moment) /i(ri) of a surviving cluster scales 
inversely with their density, /i(0) 0~^. 

The correlation length ^ close to the critical point can be determined by 
considering the linearized flow equations off criticality. ^ corresponds to the length 
scale at which the off-critical flow significantly deviates from the critical one. In this 
way, we find the famous exponential Kosterlitz-Thouless form 

ln(e/eo) - VW\ (23) 

where r = C — Cc denotes the (renornialized) distance from criticality. 

Finally, we turn to the disordered Griffiths phase (C > Cc). Because Pq flows 
to zero and Rq to infinity in this regime, we need to use the full fiow equations (jlSp 
and p9p rather than the linearized ones to determine the low-energy renormalization 
group fiow. Solving and ([TS]) in the appropriate limit gives 

Po(r) = P.exp[-(i?,/a;)e^n ^ R^{T) = e'"'^ (24) 

where Px = exp(C /x) and Rx is an integration constant. In the disordered Griffiths 
phase J <C e, the system therefore consists of essentially independent clusters with a 
distribution of characteristic energies e. Their density of states can be obtained by 
inserting (|24p into (II7p and transforming back to the original distribution p(e). This 
yields 

p(e)^exph(P,/x)(17o/e)1 (25) 

with subleading power-law corrections. We conclude that the Griffiths singularities in 
the superohmic case are exponentially weak, in contrast to the Ohmic case which leads 
to strong power-law Griffiths singularities. It is worth noting that the functional form 
of the rare region density of states (|25p within our renormalization group approach 
agrees with that predicted by heuristic optimal fiuctuation arguments [181 128j . 

5. Observables 

To find the average order parameter, we run the renormalization group to F 00 
and determine the total moment /i of the surviving clusteil^. In the disordered phase, 
C > Cc, the flow is towards Pq = and Pq = 00. Thus, all clusters are decimated in 
the late stages of the renormalization group flow. The total moment of the surviving 
clusters is therefore zero, implying that the average order parameter m vanishes, as 
expected. At criticality, the moment of a surviving cluster scales as /i(ri) ^ 
while their density behaves as n(ri) ~ fi^ (see Sec. [3]). The average order parameter 
is given by m = n(r2) /i(r2) in the limit f2 — 0. We conclude that m is nonzero at 
the critical point. To discuss how the order parameter changes upon moving into 

§ We emphasize that the strong-disorder renormaUzation group approach does not contain spin-wave 
type excitations of the ordered phase. Whether or not an infinite surviving cluster indeed supports 
long-range order depends on the strength of the corresponding transverse fiuctuations, as discussed 
in Sec. [3 
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the ordered phase, we note that the dependence oi m on r — C — stems from 
the rare renormahzation group steps that decimate clusters in the approach to the 
fixed point. By comparing the resulting losses of clusters along the critical and off- 
critical renormahzation group flows, we again find Kosterlitz-Thouless type behavior, 
m(r)— to(0) ^ ^J\r\/x. In summary, the order parameter jumps from zero to a nonzero 
value upon reaching the critical point and then increases in a square-root fashion. 

The dependence of the order parameter susceptibility on temperature T can be 
found by running the renormahzation group to — T and treating the surviving 
clusters as independent because their interactions J are small compared to T . (Note 
that this constitutes an approximation in the superohmic case because the distribution 
7r(J) does not become infinitely broad at criticality, but it becomes asymptotically 
exact for a; ^ 0.) The total susceptibility is a sum of independent Curie contributions 
from the surviving clusters, x{T) = n{T)n'^{T)/T ~ T-(i+^). 

To calculate the spin- wave stiffness ps, we apply boundary conditions at the two 
ends of the chain, x — Q and x = L, such that the rotors at the two ends are at 
a relative angle 9. In the limit of small Q and large L, the free energy density / 
depends on Q as /(6) — /(O) = ^Ps {Q/L)^ which defines ps- By optimizing the 
local twists according to the local interaction strengths, one finds ps ~ {1/J)~^ 
where (• • •) denotes the average over 7r(J). A full calculation would need to include 
the internal stiffnesses of the surviving clusters. However, we can obtain an upper 
bound for the stiffness by treating the surviving clusters as rigid. On the line of fixed 
points Pq > x,Rq = characterizing the critical point and the ordered phase, the 
distribution 7r(J) follows from transforming ([T7| back to the original variables. It 
reads tt{J) = PqU~^o ~^ . The integral dJ ti{J)/ J diverges at the lower bound 
for all Pq < \. We thus conclude that the spin- wave stiffness ps vanishes not just at 
the critical point (Pq* = x), but also in part of the ordered phase (for x < Pq < 1). 

6. Conclusions and Discussions 

In summary, we have studied the infiuence of superohmic dissipation on the quantum 
critical behavior of a disordered 0{N) (j)^ order parameter field theory. We have 
shown that the static critical behavior falls in the Kosterlitz-Thouless universality class 
independent of the exponent that controls the dissipation strength. The dynamical 
scaling is of conventional power-law type with the dynamical exponent z varying from 
1 in the dissipationless case to infinity in the limit of Ohmic dissipation. In both 
limiting cases, our theory agrees with known results. For zq = 1 we reproduce the 
dirty boson results of Ref. [12] , and in the limit zq = 2 our theory becomes identical to 
the Ohmic case [T11|T3]. In these final paragraphs, we discuss a number of remaining 
questions, and we put our results into a broader perspective. 

First, while the exponential ansatz ()17|) provides a solution of the renormahzation 
group flow equations, it is a priori not clear that the flow for a moderately disordered 
system actually follows this solution. We have therefore implemented the recursion 
relations (|S]) and (|10l) numerically and studied the renormahzation group flow for a 
broad range of parameters. A characteristic result is shown in Fig.[2l more details will 
be published elsewhere [26] . The figure shows that the fiow agrees very well with that 
predicted by our analytical solution. (Note that within the exponential ansatz ([T7|). 
the average values of C, and /3 are given by Pq~^ and Rq^ ^ respectively.) 

Second, the recursion relations (jSj) and (flUl) are valid for broad disorder 
distributions. Close to the critical fixed point, the distribution p(t) does become 
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1/(^0 

Figure 2. Numerical implementation of the strong-disorder renormalization 
group: Trajectories of the average values of C, and f) along the renormalization 
group flow. The bare interactions J are distributed uniformly between Jo and 1; 
the bare gaps are distributed uniformly between eo a-nd 1- The exponent zq = 
The dashed line is the critical trajectory predicted by I I20I1 . The right-turns of the 
lowest two curves for very small l/(x/3) are a numerical artefact. 

infinitely broad (i?o 0). The recursion relation ([TU]) for merging two clusters thus 
becomes asymptotically exact. Because essentially all renormalization steps close to 
the critical fixed point consist of merging two clusters, as discussed in Sec. SI we 
conclude that the critical fixed point is stable against subleading corrections to the 
recursion relations. To test the effects of these corrections during the early stages 
of the renormalization group, we have included them in the numerical solution. We 
found that they do not qualitatively change the flow towards the fixed points discussed 
in the present paper |26| . 

We also point out that our explicit calculations are for the large- TV limit of the 
order-parameter field theory However, the critical behavior depends only on the 
structure of the recursion relations ([5]) and (|10p . The multiplicative structure of ^ 
simply reflects second-order perturbation theory and will certainly be the same for 
finite N . The structure of (fTO)) reflects the nontrivial power-law dependence of e on 
the cluster size for superohmic dissipation. It is the same for all continuous-symmetry 
cases, as also follows from a renormalization group analysis of appropriate nonlinear 
sigma model at its strong-coupling fixed point [28] . We thus conclude that the critical 
behavior found here will be valid for all iV > 1. 

The theory developed in this paper applies to one space dimension. In 
higher dimensions, a closed analytical solution of the strong-disorder renormalization 
group appears impossible because the topology of the lattice changes under the 
renormalization procedure. However, a few qualitative aspects of the behavior in 
higher dimensions can be read off the structure of the theory. In particular, the fact 
that we obtained power-law dynamical scaling rather than activated scaling arises 
from the power-law relation ([7]) between the "gap" (distance from criticality) and the 
damping parameter of a single cluster. The damping parameter is proportional to 
the number of sites in the cluster and thus, in general, to some power of the cluster 
length scale, 7 ^ L^f where Df is the (fractal) dimension of a surviving cluster at 
criticality. Together with (O, this suggests that the dynamical exponent is generally 
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given by z — Df zq/(2 ~ zq). Interestingly, this is precisely the relation found at a 
percolation quantum phase transition [35] where Df is simply the geometric fractal 
dimension of the underlying percolating lattice. The power-law relation ^ between 
the gap and the damping parameter also implies that the Griffiths singularities are 
exponentially weak in all dimensions because, when combined with the exponential 
size distribution of rare regions, it produces an exponentially small low-energy density 
of states 18|. Even though the dynamical critical behavior in higher dimensions will 
be similar to our one-dimensional results, we do not expect the Kosterlitz-Thouless 
character to survive in higher dimensions. This follows from the observation that in 
d > 1, even the dissipationless quantum rotor system has a conventional quantum 
phase transition. 

Finally, we relate our results to the general classification of phase transitions in 
the presence of disorder based on the defect (and rare region) dimensionality [SI I18j . 
As can be seen from the calculations in Sec. 13.11 a single cluster or rare region is below 
the lower critical dimension of the problem in the case of superohmic dissipation. This 
puts the transition into class A of the classification, which implies conventional power- 
law dynamical scaling and exponentially weak Griffiths singularities. The results of 
Sec. Ulare in complete agreement with these predictions. 
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Appendix 

In this appendix, we sketch the derivation of the recursion relations ([IU|) and ([TT|) for 
the case of a strong bond. To identify and integrate out the high-energy mode of the 
two-site cluster, we diagonalize the quadratic form in its action (|9]). This yields 

Sq = [Ka\lpa{i^n)\'^ + Kbli'bi^n)]'^] ■ (A.l) 

Here, ipa and tpb are the two eigenmodes, and the corresponding eigenvalues are given 

by 



1 



Ka^b = ^[d2 + d3T \j{d2 - ds f + J| j (A.2) 

with dj — rj + \j + ^j\ujn\'^/^° . Expanding in I/J2 results in 
1 



Ka,b 



r2 + A2 + r3 + A3T ./2 + (72 +73)|c^np/"°] +0 (^-^^ . (A.3) 

As the higher eigenvalue is at least J2 above the lower eigenvalue Ka, we integrate 
out the eigenmode i/jb yielding the effective action Sq = Tj^ui^ «^a|V'a(w„)p. We now 
rescale the remaining eigenmode by defining 02 = V'a/(V'o)^^^ because we wish the 
renormalized rotor to fulfill the same large- iV constraint ((^2) = 1 as the original 
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variables. In the relevant limit J2 ^ £27 £31 the rescaling factor (-0^) approaches 2. We 
thus arrive at the renormalized action 



From the frequency term, we can immediately read off the recursion relation (jlip . 
72 = 72+73, for the damping coefficients. If we combine this recursion with the single 
cluster relation ([7]) between 7 and e, we obtain the recursion relation (ITU|) for the local 
distance from criticality. 

Alternatively (and to check the consistency of the procedure), one can calculate 
the renormalized distance from criticality directly from the two-site cluster solution 
via £2 = ''2 + A2 + r3 + A3 — J2. This requires finding the Lagrange multipliers A2 and 
A3 from the large- iV constraint equations 



We have performed this calculation in complete analogy to the Ohmic case [TS] and 
rederived the recursion relation ([TU)) for the renormalized distance from criticality 62. 
The result is identical to that obtained via the renormalized damping constant 72 
above. 
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